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Resumen 


Se evaluó el efecto de la labranza primaria sobre la conductividad hidráulica 
y porosidad de un suelo arcilloso. El experimento se realizó en el rancho “Los 
Ángeles”, municipio de Saltillo, Coahuila, México (latitud 25.1108? N y longitud 
-100.9902* O). La conductividad hidráulica in situ (Kc) fue medida con permeámetro 
de discos a tensión en la superficie (Kc,) y a 8 cm de profundidad (Kc,), mediante 
un diseño de tratamientos al azar para dos sistemas de labranza, arado de discos 
(LAD), multiarado (LM) y un testigo (NL), y cuatro tensiones de humedad: -150 (T1), 
-100 (T2), -50 (T3) y -10 (T4) mm. El movimiento del agua en perfil del suelo hasta 
el punto de saturación (Ksc) se relacionó con la porosidad. La labranza se realizó 
en junio y julio del 2009, y la Kc se midió dos meses después. Con las lecturas de 
abatimiento del agua y el tiempo se calcularon las tasas instantáneas de infiltración 
en cm h”, para Kc y Ksc. NL tuvo los valores más altos de Ksc para los dos niveles 
en que se realizaron las mediciones; mientras que en la superficie, Ksc, en LM, tuvo 
valores más bajos; LAD, tuvo el mismo comportamiento a 8 cm de profundidad 
en T4-10. Se observó también que cuando la tensión fue a -10 mm se expresa el 
funcionamiento de macroporos (> 30 um). Las velocidades para Ksc, en los cuatro 
niveles de tensión no presentó diferencia significativa, mientras que para Ksc,, 
las diferencias fueron altamente significativas (alpha = 0.05), cuyas velocidades 
promedio fueron de 87, 37 y 54 cm h* para NL, LAD y LM, respectivamente. 


Palabras clave: porosidad edáfica, permeámetro de carga variable, multiarado. 


Introducción densidad aparente (Da), distribución del tamaño 


de poros, y la conductividad hidráulica satura- 


En suelos agrícolas, la estructura de la capa 
arable es la propiedad física que más se modifi- 
ca por el laboreo mecánico, afectando con ello 
el almacenamiento y la disponibilidad de agua, 
nutrientes y oxígeno en la zona radicular, que 
es controlada por el espacio poroso del suelo 
(Buczko et al., 2006; Peth et al., 2010). Evaluaciones 
de la estructura en términos de los cambios en la 


da y no saturada son variables aceptadas como 
indicadores para valorar el efecto del laboreo 
en el estado físico del suelo (Buczko ef al., 2006; 
Peth et al., 2010); por ejemplo, Schwárzel et al. 
(2011) encontraron que la abundancia de poros 
modificados por el laboreo son inestables y 
terminan colapsados como consecuencia de los 
ciclos de humedecimiento y secado. 
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Como parte de la labranza primaria, el 
multiarado es un implemento alternativo 
al arado de discos o de vertedera, usado 
actualmente en varias zonas de importancia 
Caribe, 
debido a que está orientado a la conservación 


agrícola de Latinoamérica y el 


de suelos y del agua (Ventura et al., 2003; 
Martínez-Gamiño y Jasso-Chavarría, 2006), 
pero el efecto de este implemento sobre las 
propiedades físicas del suelo, en particular 
en su comportamiento hidrodinámico, en la 
actualidad presenta muchas interrogantes, 
pues la mayor parte de la evaluaciones 
realizadas se han concentrado en su 
desempeño mecánico (Gutiérrez ef al., 2004) 
y sobre la producción de cosechas (Martínez- 
Gamiño y Jasso-Chavarría, 2006). 

La conductividad hidráulica es una 
propiedad del suelo que indica la movilidad 
del agua y depende del grado de saturación 
y la naturaleza del mismo (Gómez-Tagle 
et al., 2008; Del Valle et al., 2009; Peth et al., 
2010), puesto que el movimiento del agua 
del suelo tanto en su estado transitorio como 
en el estacionario es un buen reflejo de la 
estructura del sistema poroso, resultante 
de la yuxtaposición de las partículas de 
suelo y de los agregados (Elrick y Reynolds, 
1992). Entre las propiedades hidráulicas 
que permiten cuantificar el tamaño, la dis- 
tribución y la continuidad del sistema po- 
roso son la conductividad hidráulica satu- 
rada, el potencial mátrico (ym) y la sortividad 
del suelo, siendo posible medir estas pro- 
piedades directo en campo mediante el uso 
del permeámetro de disco, instrumento 
ampliamente usado en la actualidad (Perroux 
y White, 1988; Regalado et al,. 2003; Moret 
y Arrúe, 2007; Gómez-Tagle et al., 2008; 
Schwaárzel et al., 2011). 

Con base en el punto anterior, el objetivo de 
la presente investigación fue evaluar el efecto 
de la labranza primaria en la calidad física de 
un suelo arcilloso, utilizando las propiedades 
hidrodinámicas, como conductividad hidráu- 
lica y porosidad. 


Materiales y métodos 


Características ecológico-geográficas del área 
de estudio 


Localización 


El experimento fue establecido en el rancho 
“Los Ángeles”, municipio de Saltillo, Coahuila, 
ubicado al noreste de México (latitud 25.1108? 
N y longitud 100.9902” O), en la subprovincia 
denominada Gran Sierra Plegada, en la parte 
norte de la Región Fisiográfica Sierra Madre 
Oriental. El clima es semiárido [BWhw(x) 
(e)], con una precipitación media anual de 350 
mm, y una temperatura media anual de 21 
“C (García, 2003). Los suelos del lugar están 
formados principalmente por suelos del tipo 
luvisol en un 40% de la superficie total (INEGI, 
1976). 


Caracterización de las parcelas experimentales 


El trabajo se realizó dentro de una sección de 35 
ha del rancho “Los Ángeles”, dedicado al uso 
agrícola desde hace poco más de 25 años, el cual 
había permanecido en descanso en los cinco 
años anteriores a la presente investigación. 
Antes de realizar los tratamientos de labranza, 
las características edáficas se determinaron en 
campo y laboratorio (cuadro 1). La vegetación 
predominante en las parcelas experimentales 
fue herbácea, compuesta principalmente por 
las siguientes familias: Laminaceas, Chenopo- 
daceas, Euphobaceas y Asteraceas. 


Medición de la densidad aparente 


Los valores de la Da inicial (Da) fueron 
determinados por gravimetría a partir de 
muestras extraídas con una barrena de núcleos 
a cuatro profundidades (5, 10, 20 y 25 cm) y la 
Da final (Da) fue calculada por el método de la 
cadena para medir rugosidad en la superficie 
del suelo, técnica descrita por Saleh (1993), 
la cual también ha sido utilizada para medir 
depresiones, erosión del suelo (Kamphorst et 
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Cuadro 1. Caracterización de las parcelas experimentales. 


Profundidad pH CE Arcilla Limo Arena Da Clase N M.O. er 
cm ds/m % e/g textura % 
0-30 8.3 ell 44 24 32 900 Arcilla 0.2 4.2 37.3 
30-60 8.3 0.51 58 28 14 1100 Arcilla 0.1 2.8 37.7 


Fuente: Laboratorio de Calidad de Agua y Suelo del Departamento de Riego y Drenaje, UAAAN; pH = potencial hidrógeno; CE = conductividad 
eléctrica en desisimens (ds/m) a 25 *C; Da = densidad aparente (kg m”); N = nitrógeno (%); M.O. materia orgánica (%); CaCO, Tot. = carbonatos 
totales en %; clase textural de acuerdo con el Sistema de Clasificación de la USDA. 


al., 2000; Merril et al., 2001) y variaciones en 
la Da por efecto de la labranza (Cadena et al., 
2004). 


Diseño experimental y análisis de datos 
Tratamientos 


El trabajo de investigación fue diseñado bajo 
un sistema completamente al azar, con tres 
repeticiones, considerando como fuentes de 
variación los siguientes tratamientos: 1) testigo 
sin laboreo mecánico, NL; 2) labranza con arado 
de discos, LAD; 3) labranza con multiarado, 
LM, y cuatro niveles de tensión: -150, -100, -50 
y -10 mm. El control de las tensiones se realizó 
de acuerdo con lo descrito por Moret y Arrúe 
(2007). 

La descripción de los tratamientos de 
labranza es la siguiente: 


1. LAD (figura la). Se realizó la labor a 0.25 
m de profundidad con un arado de tres 
discos de acero cóncavos de 60 pulgadas 
de diámetro, utilizando como fuente de 
potencia un tractor John Deere modelo 
3200 de 84 HP. El resultado de la labor es 
corte, mezcla e inversión de los estratos del 
suelo, lo que resulta en un mayor cambio 
de los agregados del suelo en el perfil. 

2. LM (figura b1). Se llevó a cabo con un 
multiarado, que es una combinación de 
varios implementos agrícolas. Presenta un 
timón tipo subsuelo, pero con la punta de 
un arado “ranchero”. Las aletas posteriores 
permiten trozar las raíces. El multiarado 


está diseñado para roturar el perfil del 
suelo en forma vertical sin invertirlo o 
mezclarlo. La profundidad de trabajo 
también se estableció a 0. 25 m. Para esta 
labor se utilizó como fuente de potencia un 
tractor New Holland modelo 6810 de 95 
HP (Martínez-Gamiño y Jasso-Chavarría, 
2006). 


El testigo fue el tratamiento de no labranza, NL. 
Instrumentación y técnica de medición 


En cada tratamiento se establecieron puntos 
de muestreo en forma aleatoria para las 
pruebas de infiltración a nivel de la superficie 
(0 cm) y a 8 cm de profundidad, utilizando 
un permeámetro de disco, el cual ha sido 
usado (Perroux y White, 1988) para medir 
propiedades hidráulicas de suelo en campo, 
permitiendo con ello identificar macroporos y 
rutas de flujo preferencial derivado de cambios 
por el majo de suelos, dispositivo utilizado en 
la presente investigación (figura 2) siguiendo el 
protocolo de campo descrito por Moret y Arrúe 
(2007). 


Análisis de datos y métodos de cálculo 


Las ecuaciones empleadas para determinar los 
flujos dinámicos y estacionarios partieron de la 
representación del modelo empírico propues- 
to por Kostiakov-Lewis (Kostiakov, 1932), el 
cual ha sido ampliamente usado en estudios 
de hidrología (Landini et al., 2007; Parhi et al., 
2007; Dashtaki et al., 2009), en particular por los 
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Figura 1. Laboreo del suelo e implementos utilizados por medio de LAD (a, la) y LM (b, 1b), y efecto de cada implemento sobre 
la rugosidad de la superficie del suelo y la vegetación (c) y (d), respectivamente. 


buenos resultados reportados (Dashtaki et al., 
2009) en suelos arcillosos, cuya expresión es la 
siguiente: 


Tex? (1) 


Donde ] es la velocidad de infiltración en 
cm h”; £, el tiempo en minutos; K representa la 
velocidad de infiltración durante el intervalo 
inicial, cuando t =1; n, el parámetro que indica 
la forma en que la velocidad de infiltración se 
reduce con el tiempo y varía de -1 a 0. 

Los parámetros del modelo se obtuvieron 
mediante el cálculo de las tasas de infiltración 
a partir del gradiente de abatimiento o lecturas 


instantáneas del reservorio con respecto al 
tiempo (01/0t), lecturas que fueron leídas 
y registradas en una plantilla cada minuto 
durante un promedio de 150 minutos para las 
18 pruebas realizadas. Los dos tipos de flujo 
del fenómeno estudiado —flujo transitorio 
(Ft) y flujo estacionario (Fe)— se identificaron 
mediante los cambios en la pendiente de las 
curvas de infiltración, [, y la lámina acumulada, 
z (figura 3). 

Una vez identificados los grupos de valores 
para Ft y Fe de cada ensayo, se separaron los 
que correspondieron al flujo estacionario 
representando Ksc, los cuales, en teoría, se 
presentan cuando la diferencia entre la tensión 
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d) Válvula 


Reservorio 


Tubo 
de 
> Mariot 


Datos técnicos del permeámetro de disco 

Diámetro interno del disco (d) = 180 mm 

Nivel del tubo de mariot y base = 7.5 mm 
Diámetro interno del tubo grande o reservorio (d) = 27.61 mm 


Fuente: Tamari y Salgado, 2004. 


Figura 2. Permeámetro de disco (a), diámetro de humedad en la superficie (b) 
y profundidad del bulbo de humedad sobre el perfil de suelo (c), y modelo del permeámetro de disco usado (d). 


Flujo Flujo 
transitorio estacionario 


Velocidad de 
infiltración (cm/h) 
Penetración (cm) 


Tiempo (horas) 


Figura 3. Modelo clásico que representa el fenómeno de infiltración e identificación de momentos de flujo (Ft y Fe). 
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en el tubo de Mariot y fuerza de succión de los 
poros es igual a cero. Los valores de Ksc fueron 
sometidos a una prueba estadística (PROC 
ANOVA) y Prueba del Rango Estudentizado 
de Tukey (HSD) en SAS, para cada una de las 
variables de respuesta estudiadas con un nivel 
de significancia de 0.05 (SAS, 2004). 


Resultados y discusión 


Cambios en la densidad aparente y porosidad 
total 


Densidad aparente (Da) 


La caracterización edáfica en cuanto a Da es 
coincidente con las características de suelo 
analizadas y clasificados con base en el 
triángulo de texturas, ya que en la literatura 
se reporta (Narro, 1994) que para la clase 
textural denominada: “Arcilla en agregados”, 
la densidad aparente se encuentra en el rango 
de 900 a 1 100 kg m”, pues de acuerdo con 
el análisis mecánico de suelo (Hidrómetro 
de Bouyoucos) de muestras tomadas en las 
parcelas experimentales, en la profundidad 
0-30 cm, las proporciones de arena, limo y arcilla 
fueron de 320, 240 y 440 g kg”, respectivamente 
(cuadro 2). 

De acuerdo con ello (cuadro 2), los resulta- 
dos indican dificultad para observar diferen- 
cias importantes en la Da, aun cuando la 


labranza es la principal fuerza que ejerce 
cambios en esta propiedad del suelo (Fuentes 
et al., 2004; Udawatta et al., 2008a y 2008b). En 
este caso se parte de una condición inicial de 
densidad muy baja, por lo que el efecto de la 
labranza en reducir la densidad es marginal, 
inclusive habría que considerar que se ha 
encontrado una relación directa entre la Da y la 
resistencia mecánica (Blanco-Sepúlveda, 2009), 
en cuyo caso, aun dentro del límite crítico 
de Da para este tipo de suelo, LM podría ser 
interpretado como un efecto de compactación 
en la capa superficial de suelo. 


Porosidad 


En complemento a lo anterior, también se calcu- 
ló la porosidad total (Pt), que teóricamente se 
puede estimar como resultado de la diferencia 
unitaria y la relación Da/densidad real, Dr 
[P+ = 1 - (Da/Dr)], donde Dr se tomó con un 
valor constante igual a 2.65 g cm” (Eynard et 
al., 2004). De acuerdo con esto, y con base en los 
valores Da antes y después de haber aplicado 
los tratamientos de labranza, se calculó la Pt, 
observando que ésta se encuentra entre 63 y 
69%, donde NL es el tratamiento que presenta 
los valores más altos (figura 4). 

Con la información mostrada (cuadro 2 y 
figura 5) se observa un efecto muy parecido 
al reportado por varios autores al evaluar el 
manejo de suelos agrícolas (Rasmussen, 1999; 


Cuadro 2. Densidad aparente (Da) y porosidad total (Pt), inicial y final por tratamiento. 


Da por estrato (cm) Pt por estrato (cm) 

Trat. 5 10 20 25 5 10 20 25 
kg m? % 

Da inicial 
NL 904 923 816 955 65.8 65.1 69.2 63.9 
LM 916 863 956 973 65.4 67.4 63.9 63.3 
LAD 973 983) 915 915 63.3 64.6 65.5 65.5 
Da final 
NL 904 925 816 955 65.8 65.1 69.2 63.9 
LM 905 852 944 961 65.8 67.8 64.4 63.7 
LAD 921 888 866 866 65.2 66.5 67.3 67.3 
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Porosidad relativa 


Profundidad en cm 


a) 


65% 70% 60% 65% 70% 
1 + 1 J 
a 
10 
T T 1 
si S 8 4 
b) 


NL AH LM —_A— LAD 


Figura 4. Perfil de porosidad por tratamiento antes (a) y después de la labranza (b). 


Lipiec et al., 2005; Cervinka et al., 2009; Lozano 
et al., 2010), puesto que la porosidad depende 
más de la organización de partículas que 
del tamaño de las mismas (Peth ef al., 2010; 
Schwaárzel et al., 2011). 

Además, se observa que Pt en LM para 
las cuatro profundidades analizadas tiene un 
incremento promedio de 1.4% + 1.9%, aunque 
en la segunda profundidad (10 cm) se presenta 
un aumento de 4.2%, mientras que LAD tiene 
un incremento promedio de 1.9% + 0.1%, 
lo cual, evidentemente, indica cambios más 
uniformes en el perfil del suelo; asimismo, se 
observa también que la labranza en ambos 
métodos (LM y LAD) tiene el mayor impacto 


en la capa superior del suelo (0-10 cm) que en 
el resto del perfil. 


Conductividad hidráulica en el suelo 
Tasa de infiltración 


Para determinar el efecto de los tratamientos 
sobre la porosidad edáfica, en los 18 sitios 
elegidos al azar se realizaron los ensayos fijan- 
do el permeámetro con tensiones equivalentes 
en mm de agua (Z1) en un gradiente de 
menos a más en el siguiente orden: -150 (11), 
-100 (12), -50 (T3) y -10 (T1), condición que 
se mantuvo hasta que cada ensayo llegó a un 
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Figura 5. Velocidad de infiltración calculada en cm h* y ajuste de polinomios de segundo grado obtenidos a partir 


de las tasas de infiltración calculadas con datos de campo medidos en la superficie (a, b y c) y a 8 cm de profundidad 
(d, e y f) para cada tratamiento de labranza (NL, LAD y LM) y cuatro tensiones en mm: -150 (11), -100 (T2), 
-50 (13) y -10 (T4). 


flujo estable. El movimiento del agua a través 
del suelo desde su fase dinámica hasta su fase 
estacionaria en términos generales se presenta 
bajo la siguiente descripción: la primera con 
sucesiones de puntos de rápido cambio, y la 
segunda con puntos con tendencia asintótica 
sobre el eje del tiempo; una vez calculadas 
las tasas instantáneas de infiltración con los 
valores del tiempo acumulado como variable 
de 
tendencia polinomiales de segundo grado, 


independiente, se obtuvieron líneas 


cuyos modelos, considerando los coeficiente 
de correlación de Pearson (R?), podrían ser, con 
excepción de LAD,, como aceptables, por ser en 
su mayoría mayores a 0.85 (figura 5). 


Modelación K-L para representar el fenómeno de 
infiltración 


Los datos obtenidos en campo fueron ana- 
lizados mediante el modelo ya descrito, 
Kostiakov-Lewis (K-L), con lo cual fue posible 
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detectar los momentos Ft y Fe, en los que se 
presentó el valor de la infiltración a saturación 
en campo Ksc en cada una de los dos niveles 
donde se realizaron las pruebas. Es así que el 
Ft en ambos niveles del suelo y las velocidades 
obtenidas fueron de 16 a 240 cm h”, y de6a 
240 cm h”, para Kc, y Kc, respectivamente; 
los tiempos promedio a los que se alcanzó la 
saturación y por tanto el Fe fue en el primero y 
el segundo fue a los 143 y 106 min (cuadro 3). 
Otras observaciones que se desprenden del 
análisis de los datos (cuadro 3) son: 1) toda 
la estructura de poros, desde los de menor 


diámetro, probablemente < 30 um (microporos) 
hasta > 30 (macroporos), sufren modificaciones 
en su morfología y, por ende, en su función 
de transporte de agua y procesos asociados, 
como intercambio gaseoso suelo-atmósfera y 
movimiento vertical de solutos; 2) los ensayos 
realizados a nivel de la superficie del suelo 
para la tensión 2 muestran tendencias muy 
semejantes para los tratamientos de labranza. 
Tanto LAD, como LM, alcanzan valores de Ksc 
entre 9 y 10 cm h”, pero los que se realizaron a 8 
cm de profundidad indican que LAD, afectó la 
porosidad, que funcionan en rangos de tensión 


Cuadro 3. Tasas de infiltración inicial y final para Ft y Fe, medidos en campo y estimados con base en 18 modelos K-L, 


con diferentes valores de R?, 


Medidos Estimados 
Trat. Iennión Modelo R? E Kc E Ksc Kc, Ksc Ft Fe 
mm K-L E E 
min cm h* min cmh? log 1 (cm h1) 
T1-150 F=120:2, 06 0.878 1 132 29 30 120 31 2.08 1.49 
T2-100 MESAS 0.6748 66 54 117 37.9 54 40 1.73 1.60 
Mt T3-50 T= 810.08 104% 0.6079 125 77 141 72 75 71 1.87 1.85 
T4-10 IIS 2 0.1089 182 180 237 134 178 143 2.24 2.16 
T1-150 T=20 516 1922 0.499 1 30 61 6.6 23 7 1.35 0.87 
T2-100 =D E 0.5406 65 18 117 84 13 10 1.10 1.01 
ae T3-50 I=1647 19768 0.4648 120 48 173 33 42 31 1.62 1.550 
T4-10 E 0.3015 77 180 211 114 134 11 2.13 2.05 
T1-150 I=13,.5879 9% 0.4252 1 24 63 2 14 4 1.14 0.63 
T2-100 T= 64.886 19:40 0.2973 64 16 128 8 12 9 1.07 0.95 
.M T3-50 1=2633.9 (02% 0.2313 127 65 187 30 55 40 1.74 1.61 
T4-10 1=35 182,200 0.2001 188 240 249 150 176 151 2.25 2.18 
T1-150 I=75.369 021 0.8117 1 9 51 32.4 75 29 1.88 1.47 
T2-100 IS WA 0.2608 52 84 95 57.6 65 56 1.81 das 
Mt T3-50 1=7796.2 4998 0.2762 101 120 134 84 116 89 2.06 1.95 
T4-10 SHIA ES 0.2052 157 240 179 174 184 162 227 22 
T1-150 T=254.06 10 0.6053 1 120 24 2 254 11 2.40 1.05 
T2-100 DETAL ae 0.2065 28 18 57 3 Y 4 0.82 0.55 
LADs T3-50 I=1643.92 1077 0.2626 61 50 97 22 73 52 1.86 1.71 
T4-10 1=19 680/40 0.8133 12% 210 183 122 176 119 2l2s 2.08 
T1-150 T=4.,5416 t01% 0.2081 1 6 35 $ 5 3 0.66 0.49 
T2-100 I= WO 0.3645 40 40 83 25 38 2 1.58 1.46 
e T3-50 1=320:12. 492 0.3649 86 9% 141 70 85 74 1.93 1.87 
T4-10 I= EDIL 0.2255 145 160 198 120 143 123 2.15 2.09 
Trat. = tratamiento; I = tasa de infiltración en cm hr; £ = tiempo; min= minutos; Kc = conductividad hidráulica inicial medida in situ; Ksc = 
conductividad hidráulica a saturación medida in situ; log =logaritmo con base 10; Kc,= conductividad hidráulica estimada; Ksc, = conductividad 


hidráulica saturada estimada. 
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equivalente a T2 y T3; y 3) los ensayos realizados 
a 8 cm de profundidad para las tensiones 
T3 y T4, con excepción de LAD,, también 
muestran tendencias de comportamiento muy 
semejantes, aunque las velocidades más altas 
son para NL, y NL,, ya que van desde 180 y 240 
cm h”, respectivamente, hasta llegar un flujo 
estable entre 134 y 174 cm h”, comportamiento 
que de manera evidente contrasta con lo 
mostrado para los tratamientos de labranza en 
ambas profundidades. 


Efecto de la labranza y la succión del suelo sobre Ksc 


Cada uno de los ensayos de infiltración inició 
poniendo el tubo de Mariot a las alturas (Z1), 
equivalentes a la succión del suelo. Para 
encontrar Ksc, considerando que durante 
un ensayo la tasa instantánea de infiltración 
inicialmente es alta y decrece en el tiempo 
hasta un valor constante llamado clásicamente 
como infiltración básica controlada por la 
conductividad hidráulica saturada (Gómez- 
Tagle et al., 2008; Del Valle et al., 2009), se 
establecieron los siguientes supuestos: 1) las 
cuatro fuerzas de succión del suelo desde T4 
hasta T1 son representativas de cuatro grandes 
grupos de poros; 2) en el perfil de suelo hay 
estratos uniformes en contenidos de humedad; 
y 3) además, existe el mismo tipo y distribución 
de partículas que conforman la capa de suelo 
estudiada. 

La interpretación de las tasas instantáneas 
de Ksc permite ver el efecto negativo de la 
labranza en la velocidad promedio para los 
ensayos realizados, principalmente a 8 cm de 
profundidad. Con respecto al tratamiento de 
referencia (NL,), las tasas de infiltración en 
LAD, en T3 y T2 fueron 20 y 4.7 veces más 
bajas que las registradas en NL; mientras que 
para LM, la diferencia fue de 1.7 y 1.4 veces, 
respectivamente (figura 6). 

Los resultados de Kung et al. (2005) 
indican que el flujo saturado (Ksc) describe 
el movimiento de agua a través de poros 
grandes (Y => 1 x 107), sin embargo, habría 
que considerar que Ksc ocurre cuando el 95% 


del espacio poroso total está lleno con agua y 
el 5% restante está lleno de aire atrapado, por 
lo cual se presume que en este caso el laboreo 
mecánico, además de los macroporos, también 
impactó microporos. 


ANOVA para Kc y Ksc, expresados como índices 
Fe/Ft 


Para detectar el impacto de la labranza sobre 
las tasa de infiltración para flujo transitorio y 
flujo estacionario, representado por las tasa 
de infiltración inicial y final, respectivamente, 
se estimaron índices Fe/Ft, por considerar 
que la labranza tiene un efecto proporcional 
sobre ambos tipo de flujo, ya que durante 
toda la prueba de infiltración, Ft es gobernado 
principalmente por los cambios que produce la 
labranza sobre los macroporos, y Fe por estar 
dominado por el funcionamiento de mesoporos 
y microporos (Y < 1 x 107), los cuales ejercen 
succiones en la matriz del suelo equivalentes 
a las tensiones del permeámetro establecidas 
entre -100 y -150 mm. 

Del análisis de varianza (ANOVA) 
resultaron diferencias altamente significa- 
tivas entre los promedios de los índices 
Fe/Ft para las estimaciones de los flujos 
(log 1), transitorio y estacionario, que se 
derivaron de las mediciones realizadas a 8 
cm de profundidad, en tanto que para las 
comparaciones entre los dos niveles (a 0 cm 
versus 8 cm), no hubo diferencia estadística 
(cuadro 4), aunque es importante mencionar 
que a este mismo nivel del suelo, aun cuando 
el efecto del laboreo estadísticamente no 
es significativo, las diferencias aritméticas 
indican que LM, tiene un mayor impacto en 
macroporos, ya que la velocidad inicial en 
T4-10 fue 1.3 veces mayor que NL, (240/180), 
y para esta misma tensión, la conductividad 
hidráulica saturada o Ksc, fue 1.1 veces mayor 
que NL, (150/134). 

A 8 cm de profundidad, en los tratamientos 
de labranza sobre la conductividad hidráulica 
se observan diferentes efectos para cada nivel 
de tensión (T4 a T1) y condiciones de humedad 
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Infiltración cm 


a] 
E 150 
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[a] 
S 
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E 100 
4 


QA 
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T1 (-150) T2 (-100) T3 (-50) T4 (-10) 2 

Tensión en mm E 

5 

Figura 6. Tasa de infiltración para cuatro tensiones en los ensayos realizados a nivel de la superficie (a) 3 

y a8 cm de profundidad (b). E 

E 

E 

en el suelo, de tal manera que de los valores 1. En T1-150, mientras que el Ff inicial en E 
o 

medidos en campo de la conductividad NL, fue de 90 cm h”, para LAD, y LM, a 
hidráulica inicial (Kc) y final (Ksc) destaca lo las velocidades fueron de 120 y 6 cm Es 
siguiente (cuadro 3): h”, respectivamente; el Fe LAD, y LM, 3 
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presentaron valores muy semejantes. En 
T2-50, mientras Ft inicial para LAD, y LM, 
fue 4.6 y 2.1 veces menor que NL,, para Fe 
en LAD, y LM, fue 19.2 y 2.3 veces menor 
que la velocidad registrada para el testigo. 
En T3-100, mientras el Ff inicial para LAD, 
y LM, son 2.4 y 1.25 menores que NL,, el Fe 
en LAD, y LM, fue 3.8 y 1.2 veces menor 
que la velocidad que alcanzó el testigo. 

2. En T4-10, mientras el Ft inicial para los 
tratamientos de labranza NE, LAD, y LM,, 
fue de 240, 120 y 160 cm h”, respectiva- 
mente, el Fe NL, presentó 174 cm h”, y para 
LAD, y LM, las velocidades presentaron 
valores muy cercanos entre ambos con 122 
y 120 cm h”, respectivamente. 

3. Los tiempos en que se presenta el flujo 
estable para las mediciones realizadas 
tanto en la superficie como a 8 cm de 
profundidad varían aproximadamente en 
la misma proporción en que se mueve Ksc, 
cuyas características son las siguientes: a) 
para la primera profundidad, NL, tiene el 
tiempo más corto con 29 min y el más largo, 
con 249 min, y corresponde a LM, para 
T1-150 para el primero, y el segundo en el 
cambio de T3 a T4, respectivamente; b) en 
la segunda profundidad, NL, se logra en 24 
min y en los siguientes cambios, de T2 a T3 
y de T3 a T4, tiene un promedio de 22 min. 


Estos resultados sobre los tiempos en que 
se logra el flujo estable tienden a estar cerca a 


lo reportado por otros autores para distintos 
tipos de suelos. La mayor parte los establecen 
en rangos de 30 a 60 min (Kosak y Ahuja, 
2005), aunque Dashtaki et al. (2009) reportan 
un máximo de 270 min; asimismo, en cuanto 
a las tasas de Ksc, Gómez-Tagle et al. (2008) 
mencionan que para un suelo arcilloso con Da 
= 1.3 g/ g, las tasas de infiltración promedio 
fueron de 223.1 cm h”, mientras que Pikul y 
Aase (2003) previamente habían reportado 
tasas en un rango de 100 a 200 cm h*; más 
recientemente, Schwárzel et al. (2011), mediante 
el uso del permeámetro de disco, reportan 
valores semejantes de conductividad hidráulica 
saturada para determinar el efecto a 7.5 cm de 
profundidad para dos años de laboreo en diez 
sitios con diferentes características edáficas. 


Conclusiones 


Aunque los pequeños cambios en la densidad 
aparente (Da) no clarificaron el efecto de la 
labranza sobre las propiedades físicas del 
suelo, es evidente que las distintas tasas 
de infiltración tanto para flujo dinámico 
como para flujo estable demostraron que el 
manejo del suelo mediante los implementos 
de labranza primaria utilizados en esta 
investigación producen cambios importantes 
en el funcionamiento del sistema poroso en la 
capa de suelo a 8 cm de profundidad, zona de 
almacenamiento del agua donde el agricultor 
pone mayor énfasis en la preparación de la 
cama de raíces para sus cultivos. 


Cuadro 4. ANOVA para los índices Fe / Ft con base en las tasas de infiltración estimadas mediante el modelo K-L de los ensayos 


realizados a nivel de la superficie y a 8 cm de profundidad. 


Tensión en mm 


Trat./prof. -150 -100 -50 -10 
NL, 0.7145 a 0.9246 a 0.9862 a 0.9618 a 
LAD, 0.6409 a 0.9244 a 0.9246 a 0.9646 a 
LM, 0.5542 a 0.8853 a 0.9228 a 0.9699 a 
NL 0.7808 a 0.9625 a 0.9454 a 0.9749 a 
LAD, 0.4381 b 0.6684 b 0.9182 b 0.9247 b 
LM, 0.7439 c 0.9191 c 0.9670 c 0.9695 c 


Letras distintas en una misma columna indican diferencia estadística (Alpha = 0.05). 
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Suponiendo una capa homogénea y to- 
mando como referencia el comportamiento 
hidrodinámico que presenta el tratamiento 
testigo, se confirma que el multiarado es un 
implemento de labranza alternativo al uso de 
arado de discos, ya que tiene menos impacto 
sobre la estructura del suelo y el movimiento 
del agua a través del perfil bajo condiciones 
saturadas. 

La técnica empleada en la presente in- 
vestigación es un buen recurso para evaluar 
cambios inducidos por el laboreo mecánico, o 
por cualquier otro medio, ya que el movimiento 
del agua en el suelo está íntimamente rela- 
cionado con las funciones primordiales de la 
porosidad edáfica; en cuanto a su densidad, 
forma y distribución en el perfil del suelo, es un 
tema para seguir incursionando en este campo 
de la investigación científica. 
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Abstract 


LÓPEZ-SANTOS, A., GONZÁLEZ-CERVANTES, G., CADENA-ZAPATA, M. és GON- 
ZÁLEZ-BARRIOS, J.L. Effect of primary tillage on the physical quality of soil, as evaluated 
by disk permeameter. Water Technology and Sciences (in Spanish). Vol. HL, No. 4, October- 
December, 2012, pp. 127-141. 


The effect of primary tillage on the hydraulic conductivity and porosity of clay soil was 
evaluated. The experiment was conducted in the “Los Angeles” ranch, Saltillo, Coahuila, 
Mexico (latitude 25.1108* N és length -100.9902* O). The hydraulic conductivity (Kc) was 
measured in situ with a tension disc permeameter at the surface (Kc) and at a depth of 8 cm 
(Kc¿). A random design treatment for two tillage systems was used, with disk plow (LAD), 
multi-plow (LM) and a control (NL), and four tension moistures: -150 (11), -100 (T2), -50 
(T3) € -10 (T4) mm. The movement of water in the soil profile up to the saturation point 
(Ksc) was associated with porosity. Tillage was conducted in June and July 2009, and Kc was 
measured two months later. Based on readings of water depletion and time, instantaneous 
infiltration rates were calculated in cm h”* for Kc and Ksc. NL had the highest values of Ksc 
for the two levels at which measurements were made; while on the surface in LM,, Ksc, had 
lower values; LAD, had the same behavior at a depth of 8 cm in T4-10. It was also noted that 
the functioning of macropores (> 30 um) is expressed when the voltage is -L0 mm. Ksc, speeds 
for the 4 voltage levels showed no significant difference, while for Ksc, the differences were 
highly significant (alpha = 0.05), with average speeds of 87, 37 and 54 cm h” for NL, LAD, 
and LM, respectively. 


Keywords: edaphic porosity, variable head permeameter, multi-plow. 
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